Connecting Direct Dark Matter Detection Experiments 
to Cosmologically Motivated Halo Models 
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Several direct detection experiments, including very recently CDMS-II, have reported signals con- 
sistent with 5 to 10 GeV dark matter that appear to be in tension with null results from XENON 
experiments. We establish the cosmologically motivated parameter space for the Galactic dark mat- 
ter velocity distribution function (VDF) and illustrate that seemingly contradictory experimental 
results with the standard analysis can be made consistent with each other within this VDF param- 
eter space. Future experimental limits should be reported for a range of well-motivated VDFs. 
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Introduction. — Dark matter comprises roughly 80% 
of the matter in the universe and seeds the formation 
of galaxies and large-scale structure. Weakly Interact- 
ing Massive Particles (WIMPs) are well-motivated candi- 
dates for dark matter, and many theoretical WIMP can- 
didates have been proposed [1-4]. Though WIMPs have 
not been detected, a variety of direct, indirect, and col- 
lider experiments are rapidly progressing in their searches 
for them [5]. 

Despite rapidly improving sensitivities and analysis 
methods, several direct detection experiments are pre- 
senting a conflicting picture. The DAMA [6], Co- 
GENT [7], and CRESST [8] collaborations have reported 
hints for low mass dark matter in the mass range ~ 
5-10 GeV. Most recently, the CDMS-II collaboration 
has reported three events in their silicon detectors that 
are not explained by known backgrounds, which when 
interpreted as a WIMP signal yields a most likely mass 
of 8.6 GeV [9]. However, all these candidate events 
arc inconsistent with the null result reported by the 
XENON100 collaboration [10]. Proposed ideas to alle- 
viate the conflict include improved characterization of 
experimental backgrounds [11, 12], particle physics ex- 
planations such as tuning the ratio of the coupling con- 
stants of WIMP scattering on neutrons and protons [13], 
or more detailed examination of the local dark matter 
density and the velocity distribution function (VDF) [14- 
16]. 

The so-called "Standard Halo Model" (SHM), which 
specifics the VDF to be a Maxwcll-Boltzmann distribu- 
tion with a cutoff at the escape velocity, is commonly 
adopted by direct detection experiments. As a conse- 
quence, uncertainties in the VDF are not a standard com- 
ponent of the analysis of experimental data. For heavy 
WIMPs, greater than ~ 20 GeV, it is relatively safe to ne- 
glect the uncertainties in the VDF because the majority 
of modern experiments are not sensitive to the variation 
of the VDF in this high-mass regime. However, for lighter 
WIMPs, which have been a large focus of recent exper- 
imental efforts, uncertainties the VDF may significantly 



affect experimental results. 

Several recent studies have shown that the Maxwell- 
Boltzmann velocity distribution is not well-motivated 
from the perspective of cosmological simulations in a 
Lambda Cold Dark Matter (LCDM) universe [17, 18], 
and different VDF models have been proposed that bet- 
ter describe simulated halos [14, 16]. Methods to present 
and to compare results from different experiments with- 
out assuming a specific VDF model have been devel- 
oped [15, 19-21], though a VDF model is still required 
to translate results from individual experiments into con- 
straints or limits on physical parameters of the dark mat- 
ter particle [22-26]. 

It has not yet become a standard in the direct detection 
community to include the uncertainties of the VDF or to 
use a VDF-independent presentation in their published 
results, possibly because the traditional vanilla WIMP 
candidate has mass of ~ 100 GeV and in this regime the 
experiments are less subject to impact of the VDF. As 
intriguing signals continue to mount, and new theoretical 
models of low mass dark matter are constructed [13, 27- 
29] , it is important to systematically address the issue of 
the VDF in the context of direct detection experiments 
and establish well-motivated models for it. In a recent 
analysis of cosmological simulations, Mao et al. [16, here- 
after Paper I] have empirically determined that the VDF 
in dark matter halos may be described by the following 
functional form with two parameters (t>o,.p), 
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(1) 
Paper I showed that a family of VDFs in the form of Eq. 1 

are consistent with current cosmological simulations, and 
also that (vo,p) depends primarily on the quantity r/r s , 
the radial position where the VDF is measured, normal- 
ized by the scale radius of the density profile of the halo. 
Due to the unknown scale radius of the Milky Way dark 
matter halo, the halo-to- halo scatter in (vo,p) from sim- 
ulations, the scatter at a fixed position within a given 



halo, and the potential impact of baryonic physics, there 
still exist significant uncertainties in the VDF, and hence 
in direct detection event rates. 

In this paper we explicitly identify the range of the 
parameter space described by Eq. 1 that is allowed by 
cosmological simulations. We further demonstrate that, 
within this well-motivated parameter domain, there is 
large range of predicted event rates for a WIMP model 
with a fiducial mass of 8.6 GeV, and show that it is pos- 
sible to simultaneously explain the CDMS-II, CoGENT, 
and XENON100 results by simply by changing from the 
SHM to a more cosmologically motivated VDF. This 
opens up the intriguing possibility that the tension be- 
tween these experiments is resolved by uncertainties in 
the Milky Way halo model, and motivates the develop- 
ment of a stronger connection between cosmological sim- 
ulations and predicted direct detection event rates. We 
conclude by discussing the current sources of the uncer- 
tainties in the VDF, the possibilities for eliminating the 
uncertainties, and suggest how to mitigate these uncer- 
tainties in the experiment analysis. 

The Distribution of the VDF Parameters. — We first 
identify a domain of allowed parameter space from cos- 
mological simulations of dark matter halos for the VDF 
of Eq. (1). Paper I identifies the best- fit VDF param- 
eters vo/v esc and p of individual halos from simulations, 
and indicates an apparent correlation between these two 
parameters for a fixed r/r s . This degeneracy between 
^o/^esc and p impedes a simple description of the param- 
eter domain of interest. To break this degeneracy, we 
instead find it useful to parameterize the VDF of Eq. (1) 
by v Tms /v e sc and p, where i> rms is the root-mean-square 
velocity, defined as 
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For notational simplicity, hereafter we use v Tms and vq to 
refer to their respective normalized values, v rms /v csc and 
v /v csc . 

In Figure 1 we show the value of i> rm s as a function 
of (vo,p)- There is an one-to-one correspondence be- 
tween (v lms ,p) and (vo,p), so the VDF of Eq. (1) can 
be completely specified by (v rms ,p). Furthermore, lines 
of constant v vms follow the relation between vo and p for 
a fixed r/r s ; v Tms is largely determined by r/r 8 , while 
the halo-to-halo scatter is primarily determined by the 
parameter p. This is physically explained by noting that 
v rms is the ratio of the average energy to the escape en- 
ergy, which is directly related to the relative position in 
the gravitational potential. 

With this parameterization, we can now specify the 
parameter domain of interest within cosmological simula- 
tions. Figure 2 shows the 90% scatter on the VDF param- 
eters for three different samples of simulated halos. One 
sample is from the Rhapsody simulation [30], in which 
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FIG. 1. Contours show the value of Urms/nesc as a function of 
{vo/v CS c,p), from the VDF model of Eq. 1. 



there are 96 halos with virial mass of ~ 10 148 M Q /i 1 . 
The other two samples are halos with virial mass of 
~ 10 14 M Q h~ 1 and of ~ 10 13 M Q /i _1 respectively, in the 
the Bolshoi simulation [31]. We use samples of halos 
with different masses in order to determine if there are 
mass trends of the VDF parameters. As shown in Paper 
I and more explicitly in Figure 2, there is no mass trend 
indicated over three orders-of-magnitude in mass, imply- 
ing that it is reasonable to apply the following analysis 
to Milky Way mass halos. 

We set the domain of interest on v Tms based on the 
current observational constraint on r/r s , which is, con- 
servatively, [0.15,1.2] [16, and references therein]. This 
then sets the domain of interest on v rms to be [0.35, 0.52]. 
Since the parameter p is not affected by r/r s , guided by 
the 90% halo-to-halo scatter from Figure 2 we set the 
domain of interest on p to be [0,3]. We note that the 
magnitude of the halo-to-halo scatter is comparable to 
the directional scatter at a fixed radius within an indi- 
vidual halo, so the above domain will not shrink even 
if one could remove the halo-to-halo scatter completely, 
given our lack of knowledge about the Earth's angular 
position. These constraints may be modified by future 
observational results or a better understanding of the 
dark matter distribution in simulated halos; we discuss 
if there are foreseeable improvements of the constraints 
on this domain of interest in the last section. For the 
sole purpose of demonstrating the impact of the VDF on 
the direct detection, one could take this full domain of 
interest to be the allowed VDF parameter space given 
modern numerical simulations. 

A Demonstration with Mock Experiments. - We 
demonstrate the impact of uncertainties in the VDF on 
direct detection experiments by considering two mock ex- 
periments, which we call Exp. X and Exp. S, and inves- 
tigate how the different parameters of the VDF in Eq. 1 
impact the interpretation of the results. In this demon- 
stration, we assume a WIMP model which has a mass 
mj m = 8.6 GeV and a WIMP-nuclcon cross section at 



0.50 



0.45 



0.40 



0.35 




0.6 


1 


" 


0.5 




- 


1 0.4 

5 






- 


3 ° 0.3 


t 




- 


0.2 






=t 


0.1 
n n 


i 



10" 



10 u 




r/r s 



r/r 



FIG. 2. From left to right, the plots show « rma /v e sc (from fitted profiles), fitted Vo, and fitted p respectively, as functions of r/r s , 
for simulated dark matter halos of three samples. The red sample consists of halos of ~ 10 1 M@ in the Bolshoi simulation, 
the green sample consists of halos of ~ 10 14 J\/q in the Bolshoi simulation, and the blue sample consists of all halos in the 
Rhapsody simulation (~ 10 1 ' Mq). The errorbars show the 90% halo-to-halo scatter of each sample. 



zero momentum transfer a = 1.9 x 10" 
spired by the recent results from the CDMS-II experi- 
ment [9]. Note that this mass and cross section are also 
consistent with the recent CoGENT analysis [32]. 

For Exp. X we consider Xenon as the target nucleus, 
with a nuclear recoil energy threshold of 6 keV, and an 
exposure of 6000 kg-days. For Exp. S we consider Silicon 
as the target nucleus, with a recoil energy threshold of 
7 keV, and an effective exposure of 7.1 kg-day, chosen 
to obtain a mean event count of 3 in the case of the 
SHM. In both experiments, to highlight the theoretical 
impact of the VDF we assume a sharp energy cutoff at 
the threshold energy, and both perfect energy response 
efficiency and energy resolution. We fix the local dark 
matter density to be po = 0.3 GeV/cm 3 , and assume 
equal WIMP coupling to the neutron and proton. We 
set the Galactic escape velocity to be 544 km/s, and take 
average speed of the Earth in the Galactic frame to be 
232 km/s. 

Given the parameters stated above, we can then cal- 
culate the predicted event rate R, which is the integral 
of the differential event rate per unit detector mass over 
the recoil energy Q, 



R = 



dQ dQ ' 



(3) 
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Here /i is the WIMP-nuclcon reduced mass, A is the 
atomic number of the nucleus, |-F(Q)| is the nuclear form 
factor [33], -u m i n = (Qmjv/2/i 2 ) 1 / 2 for an elastic collision, 
/ is the VDF in the Galactic rest frame, and v e is the 
velocity of Earth in the Galactic rest frame. 

The question we address in this demonstration is how 



the probability of a certain experiment observing TV col- 
lision events (assuming all the events are real WIMP- 
nucleus collisions) varies with different models for the 
VDF. We define Px to be the probability that Exp. X 
observes no events, and Ps the probability that Exp. S 
observes three events. We calculate the probabilities as- 
suming that WIMP-nucleon collision events follow a Pois- 
son process, 



P(N\\) = -e-\ 



(5) 



where A" is the number of events, which equals for Px 
and 3 for P s , and A is a dimensionless parameter that 
equals the predicted rate times the exposure of the ex- 
periment. Note that A changes with the WIMP model, 
the experimental setup, and the VDF. In the demonstra- 
tion we always fix the WIMP model and the settings of 
each of the two experiments, and only change the VDF 
to see its effect. 

Assuming the SHM, we obtain Px = 4.65 x 10~ 7 and 
Ps = 0.224. With these assumptions (including the 
sharp energy cutoff), given the low Px, Exp. X rejects 
the WIMP model at a high confidence level. So if Exp. S 
does indeed observe WIMP events, it implies a strong 
tension between these two experiments. Note that when 
the SHM is assumed, this conflict remains for any es- 
cape velocity larger than 515 km/s. However, the results 
change dramatically if a different VDF model is assumed. 
Assuming the VDF in Eq. (1) with a range of parame- 
ters motivated from cosmological simulations, we calcu- 
late Px and Ps and show the results in Figure 3. 

The uncertainties in the VDF can have distinct effects 
on the different experiments. Figure 3 shows that Px 
is a strong function of p, while Ps only mildly depends 
on i> rms and is insensitive to p. Because different exper- 
iments have different responses to changes in the VDF, 
a given VDF can reconcile two experiments that are in- 
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FIG. 3. The contours show the probabilities Px (left), Ps (middle), and Px x Ps (right), as functions of the VDF parameters 
t'rms/fcsc and p in the region of interests. Color scales on the three plots are all the same. Px is the probability that Exp. X 
observes no event, and Ps is the probability that Exp. S observes 3 events. Values below 0.05 are excluded with 95 percent 
confidence. High values of p can significantly reduce the tension between the two experiments, when compared to the SHM. 
See text for details of the mock experiments and other assumptions. 



consistent with one another when using the SHM. 

In the left panel of Figure 3, we see that Exp. X, which 
is strongly ruled out with the Standard Halo Model VDF, 
can only reject less than half of the parameter domain 
at a 95% confidence level when the VDF is allowed to 
vary. On the other hand, Exp. S could still observe three 
events, given that Ps > 0.05 for almost all v rms and p 
within the ranges shown on Figure 3. In the right most 
panel we also show the joint probability Px x Ps- In 
roughly one-third of the parameter domain, the possibil- 
ity of Exp. S observing three events and Exp. X observing 
none cannot be excluded. In order to exclude this WIMP 
model for all possible VDF considered within this do- 
main, we find that Exp. X must lower its energy threshold 
to at least 5.25 keV, if all other conditions and assump- 
tions unchanged. This may be possible with forthcoming 
xenon-based experiments [34]. 

As a reminder, the above analysis does not include the 
effect of background noise, the energy cutoff, the energy 
response efficiency, and the energy resolution of the mock 
experiments, and hence caution should be invoked when 
drawing strong conclusions about the regarding the re- 
lation between XENON100 and CDMS-II experiments. 
However, it does clearly motivate a full self-consistent 
statistical analysis with a VDF of the form Eq. (1) be- 
cause if the dark matter is in fact a light WIMP, a model 
for the VDF will be required to translate measurements 
into physical parameters of the dark matter particle. 

Quantifying Uncertainty from the VDF. - The re- 
sults of the above demonstration can be also interpreted 
as follows: under the assumption of the WIMP model, if 
Exp. S observes three true WIMP events and the Exp. X 
observes none, the SHM is ruled out by the experiments, 
while some range of our VDF parameter space is still al- 
lowed. Within this allowed VDF parameter space, from 
the perspective of cosmological simulations is there a sin- 



gle preferred set of the VDF parameters that could be 
adopted by experimentalists? Unfortunately, there is a 
range of uncertainties that impede a simple choice. 

Paper I presented a detailed discussion of the sources 
of scatter. Here we further distinguish these sources ac- 
cording to their contribution to the uncertainties in w rms 
or in p. We find here that v rms is largely determined by 
r/r s ; the uncertainty in this parameter is thus driven by 
observational uncertainty in r/r s for the position of the 
solar system with respect to the density profile of the 
Milky Way. Conservative estimates of the concentration 
parameter of the Milky Way imply the region of v ms used 
in Figure 3; with more optimistic assumptions one can 
constrain r/r s £ [0.32,0.50] [35]. This will narrow the pa- 
rameter range shown in Figure 3 but would not change 
our conclusions. It is likely that future data on the mo- 
tions of Milky Way halo stars and satellites will be able to 
further constrain the density profile of our Galaxy's halo 
to minimize this uncertainty. We note that the present 
analysis considers only dark matter simulations; further 
study should confirm whether the relationship between 
v T ms an d r/r s is impacted by baryonic physics. 

The uncertainty in p, on the other hand, at present 
appears to be irreducible. The halo-to-halo scatter in p 
could originate from the different intrinsic properties be- 
tween halos, but we have not yet found any significant 
correlations between p and physical properties of the halo 
(even if found, the quantity may not be well-constrained 
observationally) . In principle, one could ignore the halo- 
to-halo scatter if we had a simulation that resembles the 
Milky Way halo in every way; however, there would still 
be intra-halo scatter due to variation of VDF in vari- 
ous angular positions at a fixed radius. In Paper I, we 
found that the intra-halo directional scatter is not smaller 
than the halo-to-halo scatter. Last but not least, baryons 
could also possibly impact the shape of the VDF as char- 
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FIG. 4. The relative scatter in g(vmi n ), as denned in text, as 
a function of « m m/«csc The red solid line shows the effect 
of u lms £ [0.35, 0.52] (while fixing p = 1.5), the red dashed 
line shows the effect of a reduced parameter space w rm s £ 
[0.43, 0.46] (while fixing p — 1.5), and the blue solid line shows 
the effect of p £ [0,3] (while fixing u rms = 0.45). 



acterized by p. 

Given that there is not a single preferred set of the 
VDF parameters, when making statistical statements 
about signals or exclusions, one should marginalize over 
the VDF parameters. For example, one can calculate 
the scatter due to VDF in the differential rate, which is 
proportional to the quantity 



g{Vmin) 



,3 /(V + Ve) 
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(6) 
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where 



L(y) = min(y + v e 



, 2y, 2v e 



(8) 



and other variables are all defined in the same way as in 
Eq. (4). For the VDF of Eq. (1) and the parameter do- 
main of interest we introduced previously, we calculate 
the relative scatter of <?(v m i n ) for different v m - m . Here 
we define relative scatter to be the range (the difference 
between the maximum and the minimum) of a set of val- 
ues divided by the mean of the same set. The relative 
scatter of g(v m in) due to the uncertainty in v lms and p 
are shown separately on Figure 4. At larger w m ; n (higher 
energy bins), the scatter of <?(w m in) mainly comes from 
the uncertainty in p, and at lower i; m m, from the uncer- 
tainty in u rms , but the latter may be minimized to be 
comparable to the uncertainty in p if r/r s can be better 
determined. We note that features on the plot are due 
to the non-zero speed of the Earth in the Galactic frame. 
Since this plot shows the scatter of g(v m { n ) between halos 
or regions within the halo, these features do not show up 
on the energy spectra of the detection experiments. 

In conclusion, we show that when one assumes a single 
parameterized model for the VDF, results from different 



experiments cannot be directly compared. Even when re- 
stricting to cosmologically motivated VDFs that we have 
discussed, there is a wide range of interpretations pos- 
sible for current experimental results. Consequently, to 
present experimental results we believe the best strategy 
is to either use a VDF-independent method [15, 19, 21], 
which most conservatively only compares the region of 
Vrain that is explicitly probed by the experiments, or to 
use more different VDFs (for example, the SHM, the 
model of Eq. (1) with (vo,p) = (0.2,3), and that with 
(vq,p) = (0.215, 1.5)) to highlight the possible uncertain- 
ties. When making statistical statements about signals or 
exclusions, one should marginalize over the VDF param- 
eters; here we provide the cosmologically relevant domain 
of interest over which to do this. 
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